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Abstract. We study the perspectives of resonant charmonium production in p + A reactions within the Multiple 
Scattering Monte Carlo (MSMC) approach. We calculate the production of the resonances <I'(2S) and <? r (3770) on 
various nuclei, their propagation and decay to dileptons and D + D in the medium and vacuum, respectively, employing 
parametrizations for the D, D self energies taken from QCD sum rule studies. The elastic and inelastic interactions of 
the charmonia and open charm mesons in the medium are taken into account, too. It is found that the D, D invariant 
mass spectra from light and heavy nuclei are not sufficiently sensitive to the in-medium properties of the <I'(2S) and 
1^(3770). However, a 'suppression' of low mass dileptons from the $'(3770) might be seen experimentally as well as a 
small broadening of the <I'(2S) dilepton spectra. 
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1 Introduction 

The dynamics of charmonia and open charm mesons at finite 
baryon density and/or temperature has gained sizeable interest 
especially in the context of a phase transition from hadronic 
matter to a quark-gluon plasma (QGP), where charmed meson 
states should no longer be formed due to color screening ^ 
g]. However, the suppression of the cc vector mesons in the 
high density phase of nucleus-nucleus collisions as seen by the 
NA50 Collaboration [||| might also be attributed to inelastic 
comover scattering (cf. and Refs. therein) if the 

corresponding charmonium-hadron cross sections are in the or- 
der of a few mb [|T^,[ri||l^,[l3Hl5,|l|]. Theoretical models here 
differ in their predictions by at least an order of magnitude [ [l6| ] 
such that at present no stringent conclusions can be drawn. As 
pointed out by Seth Jl7| | this uncertainty can be substantially 
reduced by experimental data from p induced reactions on nu- 
clei, where charmonia can be formed on resonance with moder- 
ate momenta relative to the target nucleus. Such studies might 
be performed experimentally at the HESR, which is proposed 
as a future facility at GSI-Darmstadt Q. 

In previous studies we have explored the perspectives of 
charmed meson - nucleon scattering in the pd reaction [[I9|], 
where the charmonium is produced on resonance with the pro- 
ton in the deuteron and may scatter with the spectator neutron. 
We have also investigated the production of the resonances 
^(3770), 5^(4040) and !f;(4160)_on various nuclei, their propa- 
gation and decay to D, D, D* , D* , D S ,D S in the medium and 
vacuum [20|, respectively. Furthermore, the elastic and inelas- 



tic interactions of the charmonia and open charm mesons in 
the medium have been followed up in the Multiple Scattering 
Monte Carlo approach (MSMC) to study the collisional effects 
as a function of the target mass A [p0[]. In our previous stud- 
ies we have discarded the real part of the meson self energies 
dUI(M); the latter can be expressed as single-particle poten- 
tials or 'mass shifts' via dill/ (2M) and might be calculated by 
dispersion relations. 

Another way of getting information about mass shifts of 
hadrons in the medium are QCD sum rules. Related studies 
for J/W mesons suggest that the mass shift 5Ri7y/(2Afy) at 
normal nuclear matter density is only in the order of a few 
MeV | pl| , |2^ ] due to a small coupling of the c, c quarks to 
the nuclear medium. We thus will assume 3?77y « (V = 
J/<P,<I r (2S),il r (3770)) throughout the following study. How- 
ever, the open charm mesons (D, D, D*,D*) are expected to 
show sizeable mass shifts due to their light quark content [ p3[ 
^,^,^6|] in analogy to the kaons and antikaons. Especially for 
dropping D, D effective masses in the medium the decay width 
of the vector mesons then should be enhanced at finite density 
PA as suggested in Ref . [ 25 1 . 

We recall, that experiments on K ± production from nucleus- 
nucleus collisions at SIS energies of 1- 2 A- GeV have shown 
that in-medium properties of the kaons are seen in the collective 
flow pattern of K + mesons both, in-plane and out-of-plane, as 
well as in the abundancy and spectra of antikaons ||7| p7[ , ^ ,p9[ 
p(j[pl{ ]. In-medium modifications of the mesons have become a 
topic of substantial interest in the last decade trigger ed in part 
by the early suggestion of Brown and Rho []32,33|, that the 
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modifications of hadron masses should scale with the scalar 
quark condensate < qq > at finite baryon density. Thus, in 
view of the analogy to the K mesons - exchanging s, s by c, c 
quarks - in-medium modifications of the open charm mesons 
D, D should be observable, too. First exploratory studies by 
Sibirtsev et al. [ p3[ ] have suggested also an enhanced produc- 
tion of D, D mesons in p + A reactions at subthreshold ener- 
gies. Such enhanced cross sections should even more promi- 
nantly be seen in central Au + Au collisions at 20 to 25 A-GeV 
as advocated in Ref . [|4|] . However, some note of caution ap- 
pears necessary since enhanced cross sections might arise due 
to secondary reactions channels and thus do not pro- 

vide a unique signal for in-medium mass shifts. It is the aim of 
the present work to explore the possibility, if such mass shifts 
for open charm mesons might be extracted directly from dilep- 
ton or D, D invariant mass spectroscopy. 

Our work is organized as follows: In Section 2 we will 
briefly describe the input cross sections to the MSMC approach 
and evaluate the formation cross sections for the resonances 
SP(2S) and ^(3770) on protons and nuclei. The fractional de- 
cay of these resonances to open charm mesons and dileptons 
in the medium and vacuum, respectively, is calculated in Sec- 
tion 3 using 'free' and 'in-medium' spectral functions. Section 
4 concludes this paper with a summary and discussion of open 
problems. 



2 Resonance production and decay in pA 
reactions 

We here examine the possibility to measure the in-medium life 
time (or total width) of the resonances V = (&(2S), ^(3770)) 
produced in pA reactions. To this aim we describe the vector 
meson spectral function by a Breit-Wigner distribution 



A V (M 2 ) = N v 



Mr t v JM) 



(M 2 



Ml 



rni v (M)f 



(i) 

where My denotes the vacuum pole mass, Willy is the vector 
meson self energy in the medium - which vanishes in the vac- 
uum - and Ny ~ l/V (for small width r to t) is a normalization 
factor that ensures J dM 2 A(M 2 ) = 1. The total width F tot is 
separated into decay and collisional contributions as 



where 



r tot (M, PA ) = r dec {M, PA ) + r coll (M, PA ), (2) 



r dec (M, PA ) = r =( M > pa) (3) 



with r c (M, pa) denoting the partial width to the decay chan- 
nels c = DD.DD* ,D* D, e + e~,p + p,~ etc., respectively. If 
the D, D meson properties change with density pa, this will 
modify the total and partial widths accordingly. The partial 
widths for electromagnetic decays, furthermore, are assumed 
to be independent on the nuclear medium. 

The in-medium collisional width r co u (M, pa) in (Q) is de- 
termined from the imaginary part of the forward vector-meson 
nucleon scattering amplitude as 



where pa again denotes the nuclear density. Furthermore, 3?/y (0) 
is determined via the optical theorem by the total vector-meson 
nucleon cross section avN that will be specified below. 



2.1 In-medium properties of open charm mesons 

Whereas the pole masses and widths of the vector mesons My 
and open charm mesons in vacuum are rather well known [ p7| ] 
we have to model the in-medium properties that enter the spec- 
tral function ([j]) via (2) and (@). As noted before, there should 
be also shifts of the vector meson pole masses by Willy (M). 
According to QCD sum rule results [|1J the J/<P meson mass 
shift is only in the order of a few MeV at nuclear matter den- 
sity po and thus of minor importance. However, QCD sum rules 
have limited predictive power for the excited states \P(2S) and 
^(3770) and sizeable mass shifts might occur e.g. from DD 
loops in case of the ^(3770) [^Sfl. Here we do not speculate on 
such effects and will discard explicit mass shifts of the \P(2S) 
and ^(3770) and concentrate on the modification of the de- 
cay channels in the medium. To this end in each channel c the 
relative decay to mesons mi and m.2 is described by a matrix 
element (squared) and phase space, i.e. 



(5) 



4ir 

r co u = — s/V(o)p^, 

My 



(4) 



where p c is the meson momentum in the rest frame of the 
resonance. The assumption (|5|) should hold in leading order, 
however, might fail if the vector meson has an explicit internal 
structure as advocated e.g. in Ref. ||3^]. Nevertheless, in case of 
the 1^(3770) only the lowest channel to DD contributes such 
that the matrix element in (Q) can be fixed by the total width 
at half maximum in vacuum, which we take as Ifwhm = 25 
MeV p7|]. We assume the same matrix element also for the 
1 ? r (2S'), which cannot be controlled experimentally since this 
state (of mass 3.686 GeV) is below the D + D threshold. The 
resulting spectral functions S(M) — 2MA(M) - including 
the electromagnetic and light hadron decay widths - are shown 
in Fig. 1 for the W(2S) (upper part) and ^(3770) (lower part) 
for the 'free' case by the solid lines. 

In the nuclear medium the D, D mesons are expected to 
show attractive mass shifts approximately linear in the density 
Pa according to the QCD sum rule studies in Ref. [^5[. We thus 
adopt 

AM D m w -a D m{p)— (6) 
Po 

where p denotes the meson momentum relative to the nuclear 
target at rest and po w 0.168 fm~ 3 the nuclear matter density. 
The coefficients a D iq (p) should be determined by a micro- 
scopic theory, however, we will first approximate them by the 
constants an = 0.06 MeV and = 0.01 MeV, respectively, 
following the suggestions in Ref. p5[|. Since an explicit mo- 
mentum dependence of these coefficients is expected to reduce 
the self-energy effects for high momenta this particular model 
should lead to upper estimates for the medium effects. 

With the in-medium masses for D, D fixed by (^) we can 
calculate the in-medium decay widths of the charm vector mesons 
(Q) and via (Q) the total decay width (Q) assuming a total VN 
cross section of 6 mb, which is extracted from the J/H? ab- 
sorption data of the NA50 Collaboration Q]. Note that photo- 
production of J/)P's on nucleons give a total J/\P — N cross 



Ye. S. Oolubeva et al.: Interactions of charmed mesons with nucleons in the pA reaction 



3 




3.60 3.65 3.70 3.75 3.80 




3.6 3.7 3.8 3.9 4.0 



M [GeV] 

Fig. 1. The spectral functions S(M) for the !f'(25') (upper part) and 
$'(3770) (lower part) in 'free' space (solid lines) and at density po 
(dashed lines) according to the model described in the text. 

section of 3-4 mb [tttj], however, the \P(2S) and ^(3770) cross 
sections with nucleons should be larger by about a factor of 2 
due to their larger radial size. The resulting spectral functions 
S(M) — 2MA(M) - including again the electromagnetic and 
light hadron decay widths - are shown in Fig. 1 for the &(2S) 
(upper part) and ^(3770) (lower part) for nuclear matter den- 
sity po by the dashed lines. We observe substantially broadened 
spectral functions for the fy(2S) as well as ^(3770) due to col- 
lisional broadening (w 15 MeV) and a larger phase space for 
D + D decay. Note, that the threshold for D + D decay accord- 
ing to (Q) at density po moves down to 3.668 GeV - which is 
18 MeV below the &(2S) pole mass of 3.686 GeV - such that 
a larger fraction of the <P(2S) spectral distribution can decay to 
D + D in the medium. The life time of this resonance decreases 
accordingly from 656 fm/c in free space to ~ 9 fm/c at density 
po such that some fraction of the <P(2S) is expected to decay 
inside heavy nuclei such as 208 Pb. 



2.2 Production cross sections in pp and pA reactions 

The mass differential production of the vector mesons in pp re- 
actions can be described by Breit-Wigner resonance formation 
on the basis of the spectral function A(M 2 ) ([|) - employing 
the proper kinematics - provided that the branching of pp — > V 




4 5 6 7 8 

T [GeV] 



Fig. 2. The calculated cross sections for < I'(2S) and $'(3770) for pp 
(upper part) and p -f 208 Ph reactions (lower part) as a function of the 
antiproton kinetic energy T in the laboratory employing 'free' (solid 
lines) and 'in-medium' spectral functions (dashed lines) as described 
in the text. Note, that no rescattering effects are included in these pri- 
mary formation cross sections. 

is known, 

a v (M) = ^ Br pp ^ v AMA V {M) if" (A/) ( 7 ) 

where the factor 3/4 stems from the ratio of spin factors and k 
denotes the momentum of the p (or p) in the cms. Following 
our previous works JT^,p0t we adopt Br(pp — > <^(3770)) ~ 
2 x 10~ 4 which is very similar to the branching of the \P(2S) 
state . The resulting cross sections for the mesons V in pp 
reactions are displayed in Fig. 2 (upper part) as a function of 
the antiproton kinetic energy Igf= T) in the laboratory re- 
flecting the resonance structure (nl) and the kinematics from (m 
in free space. Note, that due to the small width of the !f - (25) 
(r to t ~ 0.3 MeV) the kinematics have to match the resonance 
conditions very well. 

In order to simulate events for the reaction pA — > V we 
use the Multiple Scattering Monte Carlo (MSMC) approach. 
An earlier version of this approach - denoted as Intra-Nuclear 
Cascade (INC) model - has been applied to the analysis of ri 
and oj production in pA and pA interactions in Refs. [ p4[^ ,p3l 
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[44Q . The production of the hidden charm vector resonances on 
nuclei then can be evaluated within the MSMC by propagat- 
ing the antiproton in the nucleus and calculating the elemen- 
tary pp — > V production cross section at the annihilation point, 
where the local momentum distribution is taken into account in 
the local Thomas-Fermi approximation. 

The numerical results for the formation cross sections of 
the resonances \P(2S) and ^(3770) on 208 Pb are shown in the 
lower part of Fig. 2 as a function of Tp - using the 'free' spec- 
tral functions (solid lines) - indicating that the resonance for- 
mation is further smeared out by Fermi motion and the max- 
imum in the cross section becomes much reduced in case of 
the \P(2S) inspite of the larger p annihilation cross section on 
208 Pb. As discussed in more detail in Ref. [ |20| ] the averaging 
over the Fermi motion in the target leads to kinematical con- 
ditions that do not match the resonance properties in case of a 
small total width and thus lead to a suppression of resonance 
formation. For the ^(3770) and its width of - 25 MeV the 
Fermi averaging is less crucial and compensated at Tp = 5.7 
GeV by the larger geometrical size of the Pb target. 

The formation cross section modifies only slightly for the 
tf'(3770) (dashed line) when employing the density dependent 
in-medium spectral function as described above. Thus the ex- 
citation function in the antiproton (lab.) kinetic energy Tp does 
not provide relevant information on the medium modifications 
of the 1^(3770) due to the dominance of Fermi motion in the 
nuclear target. The situation is different for the if r (2S') since 
here a large increase in the total width is expected in the medium 
(cf. Fig. 1). As a consequence the relative suppression by Fermi 
smearing becomes less pronounced and we end up with a size- 
ably enhanced cross section on nuclei for the in-medium case. 
Note, however, that \P(2S) — N dissociation reactions will re- 
duce this enhancement again as well as in-medium D + D de- 
cays such that the net effect to be seen in the e + e~ channel is 
even a ~ 40% reduction of the formation cross section relative 
to a calculation with a 'free' spectral function (see below). 



2.3 Vector meson propagation and decay 

Necessary parameters for a Monte Carlo (MC) simulation of 
rescattering are the elastic and inelastic VN scattering cross 
sections and slope parameters b for the differential elastic cross 
sections da e i / dt, which are approximated by 



of 6 mb is larger by about a factor of ~ 2 in comparison to 
the J/\P — N cross section due to the larger size of the radially 
excited states. 

The resonances are propagated with their actual mass M - 
selected by MC according to the in-medium spectral function 
([j]) - with velocity py/Ey and decay in time according to the 
differential equation 



dPy(M) 

dt 



1 



— r t v ot (M,p A )Pv(t) 



7 



(9) 



with the total width (Q), while 7 denotes the Lorentz 7-factor. 
Their decay to dileptons or DD is, furthermore, determined 
by Monte-Carlo according to the mass differential branching 
ratios at the local density pa- The off-shell propagation of the 
charmonia and evolution of their spectral functions is performed 
in line with the off-shell transport equations developed in Refs. 

The open charm mesons from the decaying vector mesons 
in the present case have a light quark q — (u, d) content apart 
from the c or c quark. In the constituent quark model we get 
D+ = (cd), D~ = (cd), D Q = (cu), D° = (cu) and the 
same composition for the related vector states. In order to es- 
timate their elastic and inelastic cross section with nucleons 
we use the analogy to KN and A'A^ interactions. It is well 
established experimentally from K + p and K~p reactions that 
light quark exchanges with nucleons ((uud) or (udd)) have dif- 
ferent strength. Whereas the K + (us) scatters only elastically 
with nucleons at low momentum - since the s cannot be ex- 
changed with a light quark of a nucleon - the K~ (su) cross 
section is dominated by resonant s-quark transfer leading to 
strange baryons such as A or 2J hyperons. Similar relations are 
expected to hold for the D,D mesons [[Z3|, where especially 
the light quark contribution should give much larger cross sec- 
tion on nucleons than the cc vector resonances. This analogy is 

18| | which, however, 



[Ha 47 



based on SU (4) fi avor symmetry 

might be broken substantially in view of the different geome- 
trical sizes of K- and D-mesons. At present this is an open 
question which has to be settled by experiment. Corresponding 
experiments have been proposed in Ref. p0|]. 

For our calculations we adopt the following cross sections 
- taken as constants in the momentum regime of interest - 



a DN 



'DN 



= 10 mb; 



d(T e l 

dt 



exp(6i), 



(8) 



where t is the momentum transfer squared. These parameters 
as well as the masses of the rescattered particles determine the 
momentum and angular distributions of the particles in the final 
state. As in our previous studies JT^,p0|] we use b= 1 GeV~ 2 
for D,D + N and VN scattering as an educated guess. 

Furthermore, the inelastic cross sections of the vector mesons 
with nucleons have to be specified. Since the relative momenta 
in the V — N system for rescattering are in the order of a few 
GeV one might use high energy geometric cross sections for 
the dissociation to open charm mesons. For the two resonances 
studied here we adopt an inelastic cross section of 6 mb as well 
as a e i = 1 mb. As mentioned above the inelastic cross section 



inel 
DN 



0; a l Sf} = 10 mb. 



(10) 



Here the inelastic cross sections of D-mesons refer to c-quark 
exchange reactions with nucleons. Furthermore, charge exchange 
reactions like D + n <-> D°p or D~p <-> D°n are incorporated 
with a constant cross section <7 qoxc . = 12 mb (cf. Ref. |[20||). 



3 Resonance production in pA reactions 

Within the model described in Section 2 we now can calculate 
the cross sections for the cc resonances in pA reactions for all 
targets of interest and obtain the information about in-medium 
resonance decays (for densities pa > 0.03/m -3 ) or vacuum 
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decays, respectively. We note that apart from the resonant pro- 
duction of open charm mesons they may also be produced di- 
rectly as (cq) and (cq) pairs in pN annihilation. This chan- 
nel strongly dominates in annihilation reactions on neutrons in 
the nucleus since the charmonium production is forbidden by 
charge conservation at low energy above threshold. At higher 
invariant energies a pion might balance the charge in the iP + tt 
final channel. For our estimates we have employed the Regge- 
model analysis of Ref. [^9[ for the elementary pN — » D, D 
cross section as also used by Sibirtsev et al. in Refs. [|23|,p4"|]. 
As already shown in Ref. [M the nonresonant production of 
D, D pairs at T/ f, = 5.7 GeV is expected to be about the same 
as the resonant production via ^(3770) excitation and decay in 
case of heavy nuclei. 

3.1 Open charm propagation and rescattering 

The open charm mesons produced by resonance decay or non- 
resonant channels rescatter in the nuclear medium either elas- 
tically or inelastically. Furthermore, they change their quasi- 
particle properties (or mass) dynamically according to (6) while 
propagating from the nuclear medium to the vacuum. This dy- 
namical problem can be well addressed by the (MSMC) ap- 
proach as well for all bombarding energies and targets. 

In order to properly interpret the final results it is instruc- 
tive to show the D + D invariant mass spectra within different 
dynamical limits. To this aim we display in Figs. 3 and 4 (up- 
per left parts) the invariant mass spectrum of D + D mesons 
from ^(3770) decays in p + 12 C and p + 208 Pb reactions at 
Tp = 5.7 GeV using the 'free' charmonium spectral function 
and branching ratios while discarding open charm and !f - (3770) 
rescattering. As seen from Fig. 3 and Fig. 4 (upper left parts) 
the in-medium decays to D + D (solid histograms) make up the 
major part of the invariant mass distribution at masses above 
3.85 GeV since the life time of high mass !f - (3770) becomes 
very short in nuclei. According to the same argument the life 
time of ^(3770) 's at low invariant masses is rather high due to 
the very limited phase space for D + D decay, such that we 
obtain a preferential decay in the vacuum (pa < 0.03 /to -3 ) 
for invariant masses below 3.8 GeV especially in case of the 
light 12 C-target. 

For the next calculations we switch on the 5^(3770) — N 
rescattering, however, discard D, D mass shifts and rescatter- 
ing (middle left parts in Figs. 3 and 4). The effects of ^(3770) 
rescatterings are esssentially a broadening of the D + D invari- 
ant mass distribution to high invariant masses and a reduction 
of the vacuum contribution (dashed histograms) due to disso- 
ciation reactions. This rather obvious medium modification is 
enhanced when including D, D rescattering (lower left parts in 
Figs. 3 and 4) which additionally leads to a depletion of the in- 
variant mass distribution in the resonance peak since rescatter- 
ing changes the invariant mass of the DD pair observed asymp- 
totically. These quite drastic medium modifications compete 
with the spectral changes of the 5^(3770) in the medium due 
to the D, D potentials. 

In the r.h.s. of Figs. 3 and 4 we show the corresponding 
D + D invariant mass distributions for the dynamical !f - (3770) 
spectral function without rescatterings (upper right parts), with 
1^(3770) rescattering (middle right parts) and also including 
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Fig. 3. The calculated D + D invariant mass distributions for p + 12 C 
reactions at Tp = 5.7 GeV within various limits for the 'in-medium' 
decay (solid histograms) and 'vacuum' decays (dashed histograms), 
respectively. The diagrams on the l.h.s. correspond to calculations 
with the 'free' spectral function for the !^(3770) whereas the dia- 
grams on the r.h.s. stem from calculations with the 'in-medium' spec- 
tral functions. In the upper plots the rescatterings of the <^(3770) 
and D, Z)-mesons have been neglected whereas in the middle plots 
!^(3770) + N dissociation has been included. The lower plots addi- 
tionally include the D, D elastic and inelastic rescattering processes 
as specified in Section 2. 



D,D interactions with nucleons (lower right parts). Whereas 
the high mass spectra differ substantially compared to the 'free' 
!f - (3770) spectral function (l.h.s.) when discarding final state 
interactions, the final mass spectra - including ^(3770) — N 
and D, D interactions - look very similar. One might have 
expected an enhancement in the low mass region below 3.77 
GeV, however, the in-medium decays to D + D at lower invari- 
ant masses finally show up in the vacuum at higher invariant 
masses (above the D + D vacuum threshold) since the open 
charm mesons regain their 'free' masses when moving out of 
the medium. 

A direct comparison of the D + D invariant mass distribu- 
tions - summing up the 'inside' and 'outside' components - is 
displayed in Fig. 5 for 12 C and 208 Pb for the case of the 'free' 
spectral function (solid histograms) and the in-medium spectral 
functions (dashed histograms). Here all final state interactions 
have been included in the calculations. Within statistics there is 
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Fig. 4. The calculated D+D invariant mass distributions forp+ Pb 
reactions at Tp = 5.7 GeV within various limits for the 'in-medium' 
decay (solid histograms) and 'vacuum' decays (dashed histograms), 
respectively. The various calculations correspond to the same limits 
as in Fig. 3. 



practically no difference between the spectra for the Pb target 
such that D + D invariant mass spectroscopy does not qualify 
much to probe the in-medium >P(3770) spectral function. This 
becomes even more apparent when comparing additionally to 
the background from nonresonant DD production (dotted lines 
'BG') which - integrated over invariant mass - is expected to 
be of comparable magnitude in cross section at Tp = 5.7 GeV 
[ 20 1 . Thus, apart from the vicinity of the resonance pole, the 
nonresonant background will dominate the high mass region of 
the invariant mass spectra. We mention that in order to calcu- 
late the invariant mass distribution from the nonresonant back- 
ground we have assumed 2-body S*-wave kinematics for the 
process pp — > DD. 

Without explicit representation we mention that the per- 
spectives to study the in-medium properties of the &(2S) by 
D + D invariant mass spectroscopy at Tp = 5.4 GeV are rather 
discouraging, since the contribution from the >T(3770) decays 
by far exceeds the signal from the tp(2S) decays (cf. Fig. 0) 
and the in-medium D, D move on-shell again, such that only 
a tiny contribution above threshold survives. According to the 
authors point of view such tiny effects will be extremely hard 
to identify experimentally. 




3.8 4.0 4.2 

M DD [GeV] 

Fig. 5. The calculated D + D invariant mass distributions for p + 12 C 
and p -f 208 Pb reactions at Tp = 5.7 GeV as obtained from summing 
up the 'inside' and 'outside' decay contributions from Figs. 3 and 4 
where all final state interactions have been included. The solid his- 
tograms correspond to the case of the 'free' <//(3770) spectral function 
while the results for the 'in-medium' spectral functions are shown in 
terms of the dashed histograms. The dotted lines (BG) show addition- 
ally the expected background from nonresonant D, D production as 
calculated by 2-body final state kinematics. 



However, it remains to be seen if the expected in-medium 
modifications of the charmonia might be detected via the dilep- 
ton decay mode since the e + e~ (or pairs do not suffer 

from strong final state interactions and also have no threshold 
in the vicinity of the charmonium pole masses. 



3.2 Resonant dilepton decays in pA reactions 

We continue with the e + e~ decay modes of the charmonia that 
have a decay width of ~ 2.2 keV for the & (2S) and ~ 0.26 
keV for the $"(3770), respectively [p7|. Inspite of these low 
decay widths the dilepton decay mode might provide essential 
information about the in-medium properties of the charmonia. 

In our actual calculations we assume the dileptonic decay 
width to be constant as a function of invariant mass M e + e - 
since in case of the charmonia the total width is very small 
compared to the pole masses. We start with !f - (3770) resonant 
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M t JGeV] M, [GeV] 



Fig. 6. The calculated e + e~ invariant mass distributions for p + 12 C 
reactions at Tp = 5.7 GeV within various limits for the 'in-medium' 
decays (solid histograms) and 'vacuum' decays (dashed histograms), 
respectively. The diagrams on the l.h.s. correspond to calculations 
with the 'free' spectral function for the ty(3770) whereas the diagrams 
on the r.h.s. stem from calculations with the 'in-medium' spectral 
functions. In the upper plots the rescatterings of the $"(3770) mesons 
have been neglected whereas in the lower plots 1^(3770) + N disso- 
ciation has been included. 



production in p + 12 C and p + 208 Pb reactions at T p = 5.7 GeV 
and consider the same limits as for the D + D decay mode. The 
numerical invariant mass spectra for these reactions are shown 
in Figs. 6 and 7, respectively. 

As seen from Figs. 6 and 7 the in-medium dilepton decays 
(40% for Pb) dominate for high invariant mass pairs whereas 
the low mass components - below the D + D threshold - only 
stem from vacuum decays (60% for Pb) in case of the 'free' 
^(3770) spectral functions (l.h.s.). This phenomenon is easily 
attributed to the long life time of the ^(3770) for low invariant 
mass such that it only decays in vacuum and 'radiates' dilep- 
tons for a long time. When including ^(3770) — TV dissocia- 
tion reactions (lower parts of Figs. 6 and 7) the shape and rel- 
ative contribution from 'inside' and 'vacuum' decays does not 
change very much, only the magnitude of the mass differential 
cross section is reduced for Pb by a factor of ^0.63. For the 'in- 
medium' ^(3770) spectral function the individual results look 
different (r.h.s.), since now the resonance can decay to DD in 
the medium at lower invariant masses, too, and the life time for 
low mass i'(3770) decreases substantially. Accordingly, there 
is also a large 'inside' fraction of dilepton decays for lower 
masses which now no longer can decay in the vacuum. The 
result is a narrowing of the i'(3770) spectral function in the 
'outside' contribution from the lower mass side. The relative 
contributions from 'inside' and 'outside' decays do not change 
substantially when including again ^(3770) + N dissociation 
reactions (lower parts of Figs. 6 and 7). Only the magnitude of 




3.6 3.7 3.S 3.6 3.7 3.8 3.9 



M M [GeV] m„ [GeV] 

Fig. 7. The calculated e + e~ invariant mass distributions forp+ 208 Pb 
reactions at Tp = 5.7 GeV within various limits for the 'in-medium' 
decay (solid histograms) and 'vacuum' decays (dashed histograms), 
respectively. The calculations correspond to the same limits as in Fig. 
6 for the 12 C target. 



the differential cross section becomes suppressed as in case of 
the 'free' spectral function (l.h.s.). 

The net spectra (possibly) to be seen experimentally are dis- 
played in Fig. 8 for both targets in case of the 'free' (solid his- 
tograms) and 'in-medium' ^(3770) spectral functions (dashed 
histograms) where the 'inside' and 'outside' components have 
been added up. Apart from the fact, that the differential cross 
sections are below 0.02 and 0.1 nb/GeV, respectively, the me- 
dium modifications only show up in a suppression of the low 
mass distributions since the high mass tails are practically the 
same. Thus the net effect is not an 'enhancement' of lower mass 
dileptons as in case of the p and u> mesons [0,^], but a relative 
'suppression'. If this suppression could be seen in experiment 
it should provide valuable information on the D, D properties 
in the nuclear medium. 

The question remains, to which extent a substantial broad- 
ening of the \P(2S) in the medium will leave its traces in the 
dilepton invariant mass spectra. Our actual calculations for the 
system p+ 12 C at Tp = 5.5 GeV show no effects within statistics 
since the average density in 12 C is substantially below nuclear 
matter density and the life time of the &(2S) is large (in fm/c) 
compared to the size of this nucleus (in fm). For p + 208 Pb at 
T p = 5.5 GeV there is a small fraction of &(2S) (about 2%) 
that decay to e + e~ inside the nucleus (dashed histogram in 
Fig. 9) and show a spectral distribution that is compatible with 
a Breit-Wigner shape of 6 MeV width (cf. thin solid line in 
Fig. 9), which is larger by a factor of 20 than the 'free' width 
of 0.3 MeV. However, when summing up the 'inside' and 'out- 
side' dilepton decays again only a small broadening in the tails 
of the invariant mass spectra survives (solid histogram in Fig. 
9). We note, that this effect might be hidden in the background 
from Drell-Yan dileptons and hard to be seen experimentally. 
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Fig. 8. The calculated e + e invariant mass distributions for p ■ 



C 



and p+ Pb reactions at Tp = 5 .7 GeV as obtained from summing up 
the 'inside' and 'outside' decay contributions from Figs. 6 and 7 where 
all final state interactions have been included. The solid histograms 
correspond to the case of the 'free' ^(3770) spectral function while 
the results for the in-medium spectral functions are shown in terms of 
the dashed histograms. 

4 Summary 

In this study we have explored the perspectives of measur- 
ing the dynamics of hidden charm vector mesons and open 
charm mesons in antiproton induced reactions on nuclei by 
D+D and e + e~ invariant mass spectroscopy. Such experimen- 
tal studies will be possible at the high-energy antiproton stor- 
age ring (HESR) proposed for the future GSI facility [jl8||. The 
Multiple Scattering Monte Carlo (MSMC) calculations have 
been based on the resonance production concept with reso- 
nance properties from the PDG [[37]] in vacuum, however, mod- 
ified in the nuclear medium according to phase space (cf. Sec- 
tion 2) involving dropping masses for the D and D mesons 

The actual values have been taken from the QCD sum rule 
studies ofRef. [g. 

We find that for a heavy nucleus like 208 Pb the fraction of 
'inside' decays to dileptons or D + D reaches ~ 40% for the 
tf'(3770), which at first sight looks promising. However, the in- 
medium effects in the D + D invariant mass spectra are hard to 




3.76 



Fig. 9. The calculated e + e invariant mass distribution for p+ 208 



Pb 



reactions at Tp = 5.5 GeV from the 'inside' (dashed histogram) and 
'outside' decay contribution (dotted histogram) employing the in- 
medium <f'(25') spectral function described in the text. The 'inside' 
component can be fitted by a Breit Wigner distribution with a width 
of 6 MeV (thin solid line). The total spectrum is given by the solid 
histogram. 



see experimentally since for invariant masses from the D + D 
threshold to the !//(3770) pole mass the final spectrum is only 
slightly modified while at high invariant masses it is covered 
by the nonresonant background (Fig. 5). The e + e+ invariant 
mass spectra from the in-medium ^(3770) show a decrease at 
small masses (below the D + D threshold) since these char- 
monia preferentially now decay via the D + D channel in the 
nucleus. Thus the net effect is not an 'enhancement' of lower 
mass dileptons as in case of the p and u> mesons [[7,50], but a 
relative 'suppression'. If this suppression could be seen in ex- 
periment it might provide valuable information on the D, D 
properties in the nuclear medium. 

Furthermore, we find that for the in-medium effects incor- 
porated the D, D decays of the &(2S) in the nucleus - that be- 
come possible due to a dropping of the D, D masses - are un- 
likely to be seen due to the much larger yield from the !//(3770) 
at the same bombarding energy (Fig. 2). On the other hand, 
a small net broadening of the <^(25) resonance tails survives 
in the dilepton channel according to our calculations (Fig. 9) 
which, however, will be hard to separate from the Drell-Yan 
background. 



We close in noting that the in-medium effects explored in 
this work are based on QCD sum rule studies for hidden charm 
and open charm mesons at vanishing relative momentum and 
low nuclear densities [ pH|25| |. In the p + A reactions investi- 
gated the mesons actually have momenta of 2-4 GeV/c [ pp| ] 
with respect to the target at rest, such that the in-medium ef- 
fects explored more likely correspond to an educated guess. 
If the D, D potentials (^) decrease substantially in magnitude 
with momentum then all effects pointed out here will decrease 
accordingly and practically no medium modifications should 
be seen experimentally. On the other hand, the open charm po- 
tentials (||) might be even larger in nature and thus the effects 
become enhanced. Nevertheless, an answer here has to come 
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from related experimental studies e.g. at the HESR [jl8f| which 
appear sufficiently promising. 
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